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Abstract 



We investigate the construction of improved actions by the Monte 
Carlo Renormalization Group method in the context of SU(2) gauge the- 
ory utilizing different decimation procedures and effective actions. We 
^ • demonstrate that the basic self-consistency requirement for correct appli- 

^ — ' cation of MCRG, i.e. that the decimated configurations are equilibrium 

£Nl ■ configurations of the adopted form of the effective action, can only be 

f*"^ ■ achieved by careful fine-tuning of the choice of decimation prescription 

CN) ' and/or action. 

O ! 

, As it is well-known, the lattice formulation is the only known non-perturbative 

' formulation of gauge theories that gives the path integral in closed form preserv- 

ing gauge invariance and positivity of the transfer matrix (unitarity). In fact, 
strictly speaking, the only known way of actually defining 'continuum' gauge 
theory non-perturbatively is by placing the lattice theory on the critical surface 
of a fixed point of infinite correlation length. Ideally, therefore, one would like to 
, have the action along the Wilsonian 'renormalized trajectory' which emanates 

from the fixed point and proceeds off the critical surface in only the relevant di- 
rections. Evolved under successive renormalization group (RG) transformations 
('block -spinnings'), this 'perfect action' could therefore be used to compute di- 
rectly at any scale (coarser lattices) without any contamination from irrelevant 
directions, and hence any regularization artifacts, and correspondingly greatly 
reduced computational effort. Concrete practical implementation of this dream, 
however, turns out to be rather difficult pQ. 

A more modest approach is based on the fact that the RG trajectory start- 
ing at any suitable lattice action will evolve asymptotically to the renormalized 
trajectory. Thus, after successive block spinnings one should in principle ar- 
rive at an 'improved action' which allows computation on coarser lattices with 
suppressed discretization errors. A way to implement this is through the use 
of the Monte Carlo RG (MCRG) in which one performs block-spinning (deci- 
mation) transformations on gauge field configurations obtained by Monte Carlo 
simulations. The basic postulate here is that the decimated configurations are 
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distributed according to the Boltzmann weight of an effective action that re- 
sulted from blocking to the coarser lattice. One, however, does not know a 
priori what this action is. Given a particular block-spinning prescription, the 
general procedure that has been followed is to assume a form of the effective 
action restricted to some subspace of possible interactions [5] , and then measure 
the couplings in this action on the decimated configurations by one of the known 
methods, the demon [3] or the 'Schwinger-Dyson' method [4]. 

The purpose of this letter is to point out that, given a choice of decimation 
prescription and a choice of an effective action, straightforward measurement 
of couplings on the decimated configurations will, in general, lead to erroneous 
results. This is because the decimated configurations will not be equilibrium 
configurations of the effective action at these couplings. Surprisingly, this basic 
requirement underlying MCRG appears not to have been enforced in its ap- 
plication to gauge theories. Careful fine-tuning of the decimation prescription 
and/or the effective action is required to satisfy this requirement, if it can be 
satisfied at all within the chosen class of decimation procedures and effective 
actions. 

We demonstrate the presence of the problem and how it can be resolved in 
SU (2) gauge theory by exploring two different decimation prescriptions: the 
Swendsen (S) decimation ,5j, and the 'Double Smeared Blocking' (DSB) deci- 
mation Both prescriptions involve a free parameter, c, which is the weight 
of staples relative to straight paths in the construction of the decimated lat- 
tice bond variable out of the undecimated lattice ones. We also explore two 
different effective action models that have been proposed in the literature: the 
multiple-representation single plaquette action [7], [8], [9] 

ON , 

plaqj = l/2 j 

and the action [TU] 

S 2 = 5>n[l - \xi/2{U P )] +X>12[1 - \xi,2{U M )], (2) 

plaq red 

containing the single plaquette (p = 1 x 1 loop) and the 1x2 planar loop in 
the fundamental representation. Here \j denotes the character of the spin-j 
representation. We use the demon method to measure couplings. 



in 


Pi/2 = Pn =2.2578 


Pi = P12 =-0.2201 


demon 
demon 


p 1/2 =2.2582(4) 
Pu =2.2576(3) 


Pi =-0.2203(3) 
P 12 =-0.2201(1) 



Table 1: Measurements of couplings by the demon method before decimation. 

To check the ability of the demon to measure couplings correctly, ensembles 
of 20000 configurations for the action (JXJ) , with jjy = 1, and ^ with couplings 
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listed in the first row of Table [T] were generated on a 8 4 lattice. The demon is 
allowed 1 sweep to set the initial energy, then 10 sweeps for each configuration 
for measurements. The results, shown in rows 2 and 3 of Table [TJ demonstrate 
that the couplings are indeed accurately reproduced, though at a greater com- 
putational cost for @. 
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Figure 1: Demon fundamental representation energy flow for DSB decimation 
and action |1| at various c values. 

Starting now with the Wilson action at coupling (3 — 2.5 on a 32 4 lattice, we 
perform decimations with scale factor of 2. Then, adopting one of the effective 
actions ((T|) or @, we proceed to measure the effective action couplings on the 
decimated configurations by unleashing the demon. We first consider action |1| 
with ]n = 5/2. Fig. [T] shows the demon fundamental representation energy as 
a function of sweeps for three values of c for a DSB decimation. The prevalent 
feature of this plot is that there is significant energy flow during microcanonical 
evolution for two of the c values shown. There is flow stabilization (equilibration) 
after about 100 sweeps. This is in fact the behavior observed for any general c 
value. Furthermore, this general flow pattern is typical for other representations 

The implication of this is clear. Suppose one, working at some chosen c value, 
measures the couplings for the effective model from the decimated configurations 
after one or a few demon sweeps (i.e. on the configurations as obtained right 
after the decimation), and proceeds to generate thcrmalizcd configurations of 
the effective action at these couplings. The decimated configurations will then 

3 




not be representative of these effective action equilibrium configurations. As 
seen in Fig. [T] the decimated configurations will evolve under microcanonical 
evolution towards equilibration at a set of different values for the couplings of 
the effective action. But by then these evolved configurations no longer are 
the true original decimated configurations obtained from the underlying finer 
lattice, and cannot be generally expected to reliably preserve the information 
encoded in the original configurations. 

Ideally, one would like to have for the measurement of couplings on the dec- 
imated configurations the same situation as that seen in the test measurement 
of couplings on the undecimated configurations (cf. Table [1] above), i.e. very 
fast demon thermalization indicating that the configurations are equilibrium 
configurations of the action for which the couplings are being measured. The 
only way out then is to seek a c, if any, for which this is realized. In the present 
case there is one such value, c € [0.065,0.067], as seen in Fig. [TJ Furthermore, 
this value shows no significant flow, indicating very rapid thermalization, also 
for the other, in particular, the adjoint representation [5]. 

A very similar state of affairs is obtained when S decimations are used in 
conjunction with |T]), but the resulting fine-tuning through selection of a c- 
parameter value is not as sharp [TT]. Overall, then, DSB is better suited for the 
action (TT]) , and can be fine-tuned so that decimated configurations are represen- 
tative of equilibrium configurations of this action. 
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Figure 2: Demon lxl and 1x2 loop energy flow for DSB decimation at 
c = 0.161. 
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Figure 3: Demon lxl loop energy flow for Swendsen decimation at various c 
values. 



We now turn to the action @. Employing DSB decimation, one again finds 
the same general picture for general values of the parameter c, i.e. significant 
flow of both the demon plaquette energy £ u and 1x2 loop energy E 12 . There 
are two special c values, c = 0.063(2) and c = 0.161(1), around which there is 
nearly no En energy flow, the second one being much more sharply defined. 
There is, however, significant flow for the E12 demon energy at both values. 
This is shown for c = 0.161(1) in Fig. [2l Similar but stronger flow for E12 is 
observed for c = 0.063(2) as well. The only c value for which E12 is constant 
is c = 0.150(2), but En shows strong flow there. It appears that there is no 
c value for which DSB decimations can be fine-tuned so that the decimated 
configurations are representative of equilibrium configurations of the effective 
action @. 

Next we consider Swendsen decimations with effective action @. Figs. Oand 
0] present typical evolution of demon En and E12 energies for various c values. 
One sees from Fig. [3] that there is approximately no En flow for c = 0.19(1) 
and c = 1.1(1). For the first value, however, there is a significant E12 energy 
evolution as seen in Fig. |4j Thus only the latter value can be used. It appears 
then that S decimation at c — 1.1 results in decimated configurations that are 
nearly equilibrium configurations of (|2j). 

For both DSB and S decimations with effective action @, the c-dependence 
of demon energies is not monotonic. E.g. for DSB decimation, En first it goes 
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Figure 4: Demon 1x2 loop energy flow for Swendsen decimation at various c 
values. 



down, then up with increasing c. Similar behavior is seen with S decimations 
(though this is not apparent in the selection of c- values shown in Fig. [3]). 
Furthermore, the direction of the variation for En is not consistent with that 
for E 12 , as can be seen from Figs. [3]and[U This is in contrast to the case of 
([1]) where the demon energies vary monotonically and consistently for various 
representations and over smaller energy ranges. Furthermore, this picture is 
stable under addition of higher representation terms. This problem with action 
([2)1 may suggest that it is in fact unstable under additions of other terms in the 
same class of interactions, e.g. other loops of length size (up to) 2 and/or other 
representations. To investigate this, and, more generally, the efficacy of ((T|) , (|2]) 



c 


couplings 


-Rlxl 


R,2x2 


R3x3 


7?4x4 


0.065 


2.5023(7),-0.3098(12) 
0.1057(16), -0.0397(16) 
0.0145(14),-0.0029(15) 


-0.0030(1) 


0.1305(9) 


0.106(3) 


-0.034(14) 


1.1 


3.2925(6),-0.2703(2) 


-0.0013(1) 


-0.0096(7) 


-0.115(2) 


-0.252(7) 



Table 2: Demon- measured couplings and difference of various size Wilson loops 
measured on decimated versus effective-action-generated configurations. First 
row: action {!]) with jn = 3 (six representations); second row: action (J2J). 



as medium to long range effective actions, we consider some observables. In 
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particular, for each of the actions (JTJ) , (J5J), wc compare N x N loops measured 
in two ways: (a) on the decimated configurations right after the decimations, 
and denoted Wff£ N ; (b) on configurations generated with the effective action 
at couplings measured after decimation at the optimal c- value (DSB for (fTJ), S 
for @), and denoted Wff™ N . The results for the difference 

= AW NxN = W£„ - Wfr$ N 
NxN w^ N W*% N 

are displayed in Tableland plotted in Fig. \5\ In contrast to action (fT]), action 
([2]) shows consistent growth of the difference Rnxn with increasing N. It is 
apparent from these results that action is failing as an accurate intermediate 
to long scale effective action, and must presumably be augmented by additional 
terms. 
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Figure 5: N x N loop differences defined in eq. ([3]). From Table [2] 



In conclusion, the basic self-consistency requirement for correct application 
of MCRG is that the decimated configurations are already equilibrium config- 
urations of the adopted form of the effective action at the couplings obtained 
from the decimated configurations. Careful fine-tuning of the decimation pre- 
scription and/or the effective action form is required in order to achieve this. 
More elaborate decimation prescriptions than the simple DSB and S may be 
defined, which involve more than one adjustable parameter, and provide more 
fine-tuning control. Furthermore, more involved actions, e.g. combination of 
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fT} and may be considered. In particular, this program is still to be carried 
out for the construction of a reliable truly improved SU (3) effective action over 
a wide range regime. 
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